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LINEAR APPROXIMATION OF THE MAX* OPERATION FOR LOG-MAP DECODING 

5 

The present invention relates to a decoder for a wireless communication 
device. 

10 Wireless communication systems are widely deployed to provide various 

types of communications such as voice and data. One such system is wideband 
code division multiple access WCDMA, whiclr has been adopted in various 
competing wireless communication standards, for example 3 rd generation 
partnership project 3GPP and 3GPP2. 

15 

To overcome data corruption that can occur during RF transmission the 
different wireless communication standards typically include some form of channel 
coding, where one common channel coding technique is turbo coding. 

20 Turbo coding involves the use of a turbo encoder for encoding a code 

segment (i.e. a data packet) and a turbo decoder for the decoding of the encoded 
code segment. 

A turbo encoder includes two convolutions encoders and an interleaver, 
25 where the interleaver shuffles (i.e. interleaves) the information bits in the packet in 
accordance with a specified interleaving scheme. 

The turbo encoder uses a first convolutions encoder to encode information 
bits within a packet to generate a first sequence of parity bits in parallel to the 
30 interleaver shuffling the information bits, where the shuffled information bits are 
encoded by a second encoder to generate a second sequence of parity bits. The 
information bits and the parity bits in the first and second sequence are then 
modulated and transmitted to a receiver. 
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5 

The information bits and the first and second sequence of parity bits are 
received by a receiver and decoded by a turbo decoder. 

The turbo decoder initially stores the received information bits and the parity 
10 bits in the first and second sequence in a buffer. Initially, the information bits and 
the first sequence of parity bits from the first convolutional encoder are retrieved 
from the buffer and decoded by a first soft input soft output SISO decoder to 
provide 'extrinsic' information indicative of adjustments in the confidence in the 
detected values for the information bits. Intermediate results that include the 
15 extrinsic information from the first SISO decoder are then stored in the buffer in an 
interleaved order matching the code interleaving used at the transmitter. 

The intermediate results, the information bits and the second sequence of 
parity bits from the second encoder are retrieved from the buffer and decoded by a 

20 second SISO decoder to provide extrinsic information indicative of further 
adjustments in the confidence in the detected values for the information bits. 
Intermediate results that comprise the extrinsic information from the second SISO 
decoder are then stored in the buffer in a deinterleaved order complementary to 
the code interleaving performed at the transmitter. The intermediate results are 

25 used in a next decoding iteration performed by the turbo decoder. The turbo 
decoder performs a predetermined number of decoding iterations before producing 
a decision on the value of the decoded information bit. 

Commonly used algorithms used within SISO decoders are the maximum a 
30 posteriori MAP decoding algorithm and the log MAP decoding algorithm. The log 
MAP decoding algorithm is analogues to the MAP decoding algorithm but 
performed in the logarithmic domain. 
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5 The MAP decoding algorithm uses forward state metrics, commonly referred 

to as alphas a, and backward state metrics, commonly referred to as betas p, to 
determine soft output results, where the forward state metrics a and backward 
state metrics p characterise a state in a trellis structure. 

10 The MAX* function is used within the log-MAP algorithm and is represented 

by MAX*(a(n),b(n)), where a(n) and b(n) are inputs to the MAX* function. The 
inputs a(n) and b(n) can be forward state metrics, backward state metrics or a 
combination of both. 

15 The MAX*(a(n),b(n)) function is equal to MAX(a(n),b(n)) plus a correction 

value where the correction value is equal to log(1+exp(-|a(«)-ft( w )|))- 

The MAX(a(n),b(n)) term of the equation is usually straight forward to 
calculate, however the correction value is relatively complicated to calculate and is 
20 usually approximated using either a linear approximation, a step approximation or 
a look-up table. 

As the state metric calculations are performed within the SISO decoder the 
values within the accumulated path metrics can overflow leading to incorrect 
25 results. 

One solution to the overflow problem involves the use of modulo arithmetic. A 
modulo n operation on a number provides the remainder when the number is 
divided by n, for example 10(binary 1010)modulo 8 = 2(binary 010) and 28(binary 
30 11100)modulo 16 = 12(binary 1100). Consequently, as can be seen from the 
examples, to determine a value for a modulo operation where the remainder is a 
value to the power of two is simply a question of masking off any unwanted bits. 



WO 2005/055435 



4 



PCT7IB2004/004420 



5 The modulo function, as illustrated in figure 1, can be regarded as a sawtooth 

function. 



An alternative implementation of the modulo function can be defined by: 



xmodF = x-2F 



, which allows negative numbers to be accommodated. 



10 This function is illustrated in figure 2. 

It is desirable to have an apparatus and method for generating a linearly 
approximated MAX* log MAP algorithm that operates on modulo functions. 

15 In accordance with a first aspect of the present invention there is provided a 

decoder for a wireless communication device according to claim 1 . 



In accordance with a second aspect of the present invention there is provided 
a method for generating a MAX* value according to claim 8. 



20 



An embodiment of the invention will now be described, by way of example, 
with reference to the drawings, of which: 

Figure 1 illustrates a graphical representation of a first modulo function; 

25 

Figure 2 illustrates a graphical representation of a second modulo function; 



30 



Figure 3 illustrates a graphical representation of the variation in the MAX* 
correction term versus \a(n)-b(ri)\ ; 
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5 Figure 4 illustrates a decoder according to an embodiment of the present 

invention. 

The curve A in figure 3 illustrates the correction term for the MAX*(a(n),b(n)) 
function (i.e. MAX*(a(n),b(n)-MAX(a(n),b(n)) as a function of \a(ri)~b(n)\), where 

10 \a(n)-b(ri)\ is the absolute value of the difference between a(n) and b(n)). 

As can be seen from curve A the correction term is greatest for low values of 
\a(ri)-b(n)\ and gradually decreases to zero as \a(ri)-b(ri)\ increases. 

15 As stated above, an easy technique for approximating the correction term is 

the use of linear approximation, as illustrated by line B in figure 3. As illustrated, 
the linear approximation provides a close approximation for the correction term for 
low values of\a(n)-b(n)\. The intersection of the line B on the \a{n)-b(ji)\ axis 

indicates the \a(n)-b(n)\ value above which the linear approximation correction 

20 term goes to zero. Consequently, using linear approximation, the intersection point 
determines a threshold value, designated C, for determining if a correction value is 
to be applied to |a(w)-Z>(w)|, where the intersection point is defined by the linear 
approximation equation. 

25 The use of the linear approximation technique allows easy calculation of the 

MAX* function, as described below. 

One suitable linear approximation equation (i.e. the correction term used) is 
given by MAX(0, H n h b ( n })/y 

30 

Consequently, the MAX* function can be written as: 
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MAX(a(n),b(n)) + MAX(0, ( c M»)-*b%) 
=MAX(0+MAX(a(n),b(n)) f ( C ~ K")~ 6 W)/ +MAX(a(n),b(n)) 
=MAX(a(n),b(n), ( C ~ H n h b ( n }J/ +MAX(a(n),b(n)) 



=MAX(a(n),b(n), + + 



To minimise the problem of accumulated state metric overflow, as discussed 
above, the above terms are converted into their corresponding 'modF' values 
15 where F is selected such that |a(«)-£>(n)| <F. F is chosen by analysing the 
algorithm and determining what would be the maximum possible value of 
|a(»)-Z>(w)j for any a(n) and b(n) that can enter the MAX* function. 



10 



To ease the hardware implementation for handling the modulo value F is 
20 preferable a value to the power of two. 

The modF of a(n) becomes a(n)modF. 

The modF of b(n) becomes b(n)modF. 

25 

However, the modF of (a(n)+b(n)+C)/2 is 



(a(n)modr i (K")modF-a(n)modF)modF + C 



modF 

J 
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5 and not ^)^+W-*i r +C) mBdjr . 

This is demonstrated by the following: 

The equation xmodF is equivalent to x-2F 
10 the floor of 'x\ 

Accordingly: 



( ^ZF^ ) ' where the L X J term is 



(a(n)+C)modF = (a(n)modF+C)modF 
15 and 

(a(n)-b(n))modF = a(n) - b(n) if and only if \a(n)-b(n)\ <F 
Using these two identities proves: 

= (a(„)m<*F + W") mMF ~ a ( n ^ F ^ F+C y)moiF. 

25 
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5 For values of C<F/2 an alternative implementation of the modulo of a linear 
approximation of a MAX function is equal to: 

^(^)mod^ + C)xnod^ + ^)modFj modF + p . 

where s is calculated from the binary expression s=[a(m) xor b(m)] and [((a(m) 
10 xor a(m-1)) and ((b(m) xor b(m-1)], where a and b are represented by m bits so 
that a(m) is the most significant bit of a and a(m-1) is next to the most 
significant bit 

This algorithm is easy to calculate in silicon as s involves only binary 
15 operations and F is chosen to be a power of two. 

A decoder 400 for implementing the above MAX* equation is shown in 
figure 4 and is arranged to output MAX(a(n)modF,b(n)modF) when 
\a(n)-b(n)\'\$ greater than the threshold value C and to output 

20 (a(n)modF + 

less than the threshold value C. If \a(n)-b(n)\ equals C then either 
MAX(a(n)modF,b(n)modF) or 

( f\ ((b(n)modF-a(n)modF)modF + C)\ Jn , 4 

I a{n)modF+ " v ' ^JmodFcan be output. 



25 The decoder 400 Includes a first substracting unit 401, a second 

substracting unit 402, a calculator 403 in the form of an adder unit and a 
selector 404 in the form of a multiplexer unit. The first substracting unit 401 , the 
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5 second substracting unit 402 and the adder unit 403 are each arranged to 
receive a(n)modF, b(n)modF and the threshold value C. 



The first substracting unit 401 is arranged to generate the sign of 
(b(n)modF-a(n)modF-C)modF. The second substracting unit 402 is arranged to 
10 generate the sign of (a(n)modF-b(n)modF-C)modF. The adder unit 403 is 

arranged to generate (afr)mod F + g&Vgdf - a{n)modF)modF + C)j modF ^ 
where the division by two corresponds to a shift in bit position by one. 



The modF operation is performed by ignoring the overflow (i.e. the carry bit 
15 of the msb bit addition is ignored). 



The output from the first substracting unit 401 , the second substracting unit 
402 and the adding unit 403 (i.e. the sign of (b(n)modF-a(n)modF-C)modF, the 
sign of (a(n)modF-b(n)modF-C)modF and 

20 (*(»)modF + «*»)mod J- -«(»)mod^)modF ± C) ^^ respective|y) 

provided to the multiplex unit 404 along with the values of a(n)modF and 
b(n)modF. 



The multiplexer 404 is arranged to output a MAX*(a(n)modF, b(n)modF) 
25 equal to a(n)modF when the sign of (a(n)modF-b(n)modF-C)modF is positive 
and the sign of (b(n)modF-a(n)modF-C)modF is negative. 
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The multiplexer 404 is arranged to output a MAX*(a(n)modF, b(n)modF) 
equal to b(n)modF when the sign of (a(n)modF-b(n)modF-C)modF is negative 
and the sign of (b(n)modF-a(n)modF-C)modF is positive. 



The multiplexer 404 is arranged to output a MAX*(a(n)modF, b(n)modF) 
10 equal to ^n)mcdF j M ^ F -^ F ^ F+C) ) ^^ *e 

sign of (a(n)modF-b(n)modF-C)modF is negative and the sign of (b(n)modF- 
a(n)modF-C)modF is negative. 



It will be apparent to those skilled in the art that the disclosed subject matter 
15 may be modified in numerous ways and may assume many embodiments other 
than the preferred forms specifically set out as described above, for example the 
above embodiments could be arranged such that the modulo for other linear 
approximation equations can be calculated and an additional substracting unit 
could be used to determine the sign of a(n)-b(n) to assist the selection process. 



